16.451 Lecture 24: The Nuclear Shell Model Dec. l, 2005 


If we choose the right potential function V N (r), then the wave function for the whole 
nucleus can be written as a product of the single particle wave functions for all A 
nucleons, or at least schematically: 



oversimplification here... actually, it has to be written as an 
antisymmetrized product wave-function since the nucleons are 
identical Fermions - the procedure is well-documented in 
advanced textbooks in any case! 

With total angular momentum given by: 



And parity: 



Always + for an even 
number of nucleons... 



What to use for V N (r)? - three candidate potential functions: 
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Square well 

r<R: V{r)=-V Q 
r>R : V'(/')=0 


Advantage: easy to write down 
Disadvantages: 

numerical solutions only 
edges unrealistically sharp 


Harmonicoscillator 


V(r) = \Mu) 2 r 2 


Advantage: easy to write down and 
can be solved analytically 
Disadvantage: potential should 
not go to infinity, have to cut off 
the function at some f inite r and 
adjust parameters to fit data. 


Saxort-Woods 



Advantage: same shape as measured 
charge density distributions of 
nuclei, smooth edge makes sense 

Disadvantage: 

numerical solution needed 



Comparison: Harmonic Oscillator versus Woods-Saxon solutions: 3 


• since both potentials are spherically 
symmetric, the only difference is in 
the radial dependence of the wave 
functions 


• amazingly, when parameters are 
adjusted to make the average 
potential the same, as shown in the 
top panel, there is remarkably 
little difference in the radial 
probability densities for these 
two potential energy functions! 


• this being the case, the simplicity 
of the harmonic oscillator potential 
means that it is strongly preferred 
as a model for nuclei 



Figure 2-Z2 The square of the wave function limes r' for the harmonic: oscillator and 
the Woods-Saxon potential are plotted in units of fm“ *. 



Evidence that this works: 
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electric charge density, measured via electron scattering: 


z 

P(r) = e^|^,(r )| 2 

Z= 1 


charge density difference 
between 205 TI and 206 Pb is 
proportional to the square 
of the wave function for 
the extra proton in 206 Pb, 
i.e. we can actually measure 
the square of the wave 
function for a single proton 
in a complex nucleus this 
way! 



Various potential shapes lead to similar patterns of energy gaps, e.g.: 
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But the magic numbers are wrong © ! 


N/Z = 2, 8, 20, 28, 50, 82, 126 


Something else is needed to explain 
the observed behaviour... 


Reasonable nuclear 



Solution: the "spin-orbit" potential 
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• Meyer and Jensen, 1949: enormous breakthrough at the time because it was the 
only explanation for the observed pattern of "magic numbers" and paved the way 
for a "periodic table" of nuclei ... and the Nobel prize in physics, 1963! 


• simple idea: 


(http://www. nobel.se/physics/laureates/1963/index. h tmi 

- see Maria Goeppert-Meyer's Nobei Lecture link on this page) 



expectation value 


as in the calculation of magnetic moments, we can write: 


(^•s) - L(j 2 -e 2 -S 2 ) - y(j(j+l) -£(£+1) - s(s+l)) 


• but there are only two ways the orbital and spin angular momentum can add for a 
single particle nucleon state: 

a) "stretched state" j = / + j : = = y 


b) "jack-knif e state" j = / - 



s(7 + 1) 


-(1 + 1 ) 
2 



Spin-orbit force, continued: 7 


• The energy shift due to the spin-orbit interaction is between states of the same 
/ but different j; 



• the splitting is proportional to / and so it increases as the energy 'mcreases for the 

single particle solutions to V(r) 

• each state can accommodate (2j+l) neutrons or protons, each with different mj 

• empirically, the sign of the spin-orbit term for nuclei is opposite to that for atoms 

and the effect is much stronger in nuclei - the phenomenon has nothing to do with 
magnetism, which is the origin of this effect in atoms, but rather it ref lects a basic 
feature of the strong nuclear force. 

• with these features, the spin-orbit potential is the "missing link" required to 

correctly predict the observed sequence of magic numbers in nuclear physics © 
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Iff! 
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Reasonable nuclear 
plus spin-orbit. 



Accumulated 

occupancy 




Applications of the Shell Model 


Generic pattern of single particle states 
solved in a Woods-Saxon (rounded square 
well) potential model with appropriate 
spin-orbit interaction to reproduce the 
observed "magic number" pattern: 

State labels: H £ j 

where n labels the order of occurrence of 
a given / value, and the state labels for 
orbital angular momentum are: 
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Each state can hold (2j +1) neutrons 
and (2j+l) protons, corresponding to 

2(2j+l) distinct conf igurations of 
identical nucleons (m + , mj) to satisfy 
the Pauli exclusion principle 
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Quantum numbers for a nucleus: 
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First of all, consider a "closed shell", which corresponds to a completely filled 
single-particle state, e.g. ls 1/2 , lp 3/2 , etc... containing (2j+l) protons or neutrons: 

The total angular momentum is: 


l 

II 

M 

^ . i 
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V* 

M = rrij - 0 
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m 

all have the same 

each m value is different, 

j in a given shell 

ranging f rom -j to +j 



There is effectively only one conf iguration here, with total z-projection M = 0. 
Therefore, the total angular momentum of a closed shell must be J = 0! 


The total parity is: 


- - r 

— 

i (-«' - 

(- 1 / 

2j + l 

= + 


Since (2j+l) is always even, the parity of a closed shell is always positive. 



Valence Nucleons: 


11 


For a closed shell + n nucleons, the angular momentum and parity is determined by 
the n "valence" nucleons, since the closed shell contributes J n = 0 + : 


n 



<-u 

II 

M 

^4 

II 

(~i) £ 

n _ 

i= 1 




The parity is uniquely determined, but there may be several different values of J 
that are consistent with angular momentum coupling rules. Residual interactions 
between the valence nucleons in principle determine which of the allowed J has the 
lowest energy - we can't predict this a priori but can learn from experiment. 


"Holes" — for a state that is almost full, it is simpler to consider angular momentum 
coupling for the missing nucleons than for the ones that are present: 


result for 
a closed shell 





rt 

2j+l 

04 
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i=(n+l) 


magnitudes of the two partial sums have 
to be the same, and M values opposite. 


Total angular momentum for a collection of missing particles, i.e. holes, is the same 
as for that same collection of particles in a given shell model state. 



Magnetic Moments: 


12 


As in lecture 20, we can write: 


n = g } J jU N = (ju»J 


(J+l) 


Predictions are only simple in the limit of one valence nucleon, and then we have: 


valence proton: 


valence neutron: 




= Us,. 


s + 




Mn 


M 9s,n S Ml V 


Simplifying the results, we find: 



with 9 sp "h 5.58, 9 sn 3.83, 9 Ip 9 in ^ 



Application with one valence nucleon: 17 0 
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Excited states of 17 0 can be understood by promoting the valence neutron: 
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First excited state: 
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ground state full to here 



Excited states of 17 0: 
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Next excited state: J n - \ - 

-> explained by promoting a neutron from 
the filled 1 p 1/2 level to the ld 5/2 level 




Excited states of 17 O 
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Two valence nucleons: (enough already!) 
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This problem is much more complicated! The inner "core" nucleons couple to J 71 = 0 + 
but in general there is more than one possibility for the angular momentum coupling 
of the valence particles. 

a) (pp) or (nn) case: Z and N are both even in these cases, so we know that 

the ground state conf iguration will be 0 + no matter what 
shell model state they are in. Excited states will have higher 
angular momentum, with possibilities restricted by the 
Pauli principle. 

b) (np) case: Z and N are both odd in this case. Only 6 examples in the whole 

nuclear chart!!! In isolation, (np) prefers to form a bound 
state - the deuteron - with = 1 + 



Example: 18 F, Z = 9, N = 9 

np pair: no restrictions on total J. 
possibilities are 0, 1, 2, 3, 4, 5 

Ground state of 18 F is 1+ i.e. deuteron 
quantum numbers! 

N.B. valence nucleons interact with each other, or 
all J values would be degenerate! 



Magnetic moments revisited: 18 

Single particle predictions: 

j=e + l/2 n = [ g, (j — 1/2) + g s l 2 ] JU N 

j = e-i/2 ft = \g e j ( j + 3/2)/(j + i) - g s j / 2(j +i)\fi N 




Odd neutron nuclei: 19 



What is wrong? 

• the single particle model is too simple - nucleons interact with each other 

• conf igurations may be mixed, i.e. linear combinations of different shell model states 

• magnetic moments of bound nucleons may not be the same as those of free nucleons... 



Title page of a research monograph, Oxford, 1990: 


From the preface: 


Modern Theories of 
Nuclear Moments 


B. CASTEL 

Queen's University , Kingston , Canada 
AND 

I. S. TOWNER 

Chalk River Nuclear Laboratories , Chalk River , Canada 
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Nuclear moments have had a long history which parallels the develop- 
ment of nuclear physics as a whole. In 1937, Schmidt proposed a standard 
for magnetic moments based on the assumption that the free-nucleon 
magnetism persists in heavy nuclei. Then in the early 1950s attempts were 
made to explain deviations from the Schmidt estimates by invoking 
meson -exchange corrections or by modifying single-particle wavefunc- 
tions so as to include many-body effects. The early 1970s saw the 
emergence of fundamental problems in trying to reconcile quadrupole 
transitions and moments in light nuclei with shell- model predictions. 
Thus the study of nuclear moments could prove an excellent pedagogical 
tool to acquaint oneself with the complexities and elegance of some of 
the most current and powerful nuclear models. 



Even the most sophisticated nuclear models are not completely successful! 
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Fig. 2.1 Nuclear moments in the odd indium isotopes measured with laser 
spectroscopy techniques: (1) quadrupole moments of 1 = 9/2 + ground states; (2) 
magnetic moments of 7 = 9/2 + ground states; and (3) magnetic moments of 
7=1/2“ isomeric states. The broken lines are theoretical calculations based on 
core-particle coupling models (see Chapters 5 and 7) using effective coupling 
constants given in the figure. The Schmidt limit is shown as a dotted line. 




An interesting puzzle from beta decay: 
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last homework assignment: n C —> n B + e + + v e log(Ft) = 3.6 1 i = 20.4 min 

14 C _> + e - + v log(Ft) = 9.0 U = 5370 yr 

why so different? Both decays should have M if = 1 in the shell model (superal lowed) 
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decay is fast, as expected 


decay is anomalously slow. - 


permits the use for 14 C dating! 


Explanation lies in interactions between nucleons that are not included in the shell model 
but required to fit experimental details. Even so, 14 C is not satisfactorily understood. 



